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Abstract 

We formulate an existence theorem that states that given locahzed scalar and vector time- 
dependent sources satisfying the continuity equation, there exist two retarded fields that satisfy 
a set of four field equations. If the theorem is applied to the usual electromagnetic charge and 
current densities, the retarded fields are identified with the electric and magnetic fields and the 
associated field equations with Maxwell's equations. This application of the theorem suggests that 
charge conservation can be considered to be the fundamental assumption underlying Maxwell's 
equations. 
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I. INTRODUCTION 



The traditional presentation of Maxwell's equations follows the historical approach in 
which electrostatics and magnetostatics are studied first. Then Faraday's induction law 
is introduced to consider quasistatic phenomena. Finally, we introduce the displacement 
current to insure charge conservation and obtain the Ampere-Maxwell law, which completes 
the set of equations known as Maxwell's equations. 

The reader might wonder why an axiomatic approach to Maxwell's equations is not usu- 
ally presented in undergraduate textbooks, although an axiomatic presentation of quantum 
mechanics and general relativity, for example, can be found. ^'^ One virtue of the axiomatic 
approach is that it offers the shortest way to the essence of a theory and enables a more 
rigorous formulation.^ The basic problem for an axiomatic presentation of Maxwell's equa- 
tions is recognizing the essential postulates underlying Maxwell's equations. One of these 
postulates is unavoidable: charge conservation, which is expressed by the continuity equa- 
tion for the charge and current densities. As is well-known. Maxwell's equations imply the 
continuity equation, but is the converse implication true? If the continuity equation implies 
Maxwell's equations, then charge conservation should be considered as the fundamental 
axiom underlying these equations. 

In this paper we show how Maxwell's equations can be obtained from the continuity equa- 
tion. We formulate an existence theorem that states that given localized time-dependent 
scalar and vector sources satisfying the continuity equation, there exist two retarded fields 
that satisfy a set of four coupled field equations. When the theorem is apphed to the usual 
electromagnetic charge and current densities, the retarded fields are identified with the elec- 
tric and magnetic fields in the form given by Jefimenko^ and the associated field equations 
are naturally identified with Maxwell's equations. Therefore, not only do Maxwell's equa- 
tions imply the continuity equation, but the continuity equation implies Maxwell's equations. 
The application of the theorem to electromagnetic sources suggests an axiomatic approach 
to Maxwell's equations in which charge conservation is considered to be the fundamental 
axiom underlying these equations. 
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II. EXISTENCE THEOREM FOR TWO RETARDED FIELDS 



In this section we formulate and demonstrate the following existence theorem: Given the 
localized sources p(x, and J(x, which satisfy the continuity equation, 



there exist retarded fields F(x, t) and G(x, t) defined by 



that satisfy the field equations 



dp 

or 

'at 



Rc^ 



'at 



(1) 



(2a) 
(2b) 



V • F = ctp, 

V • G = 0, 



V X G 



a dt 



(33. 



(3a) 
(3b) 

(3c) 
(3d) 



The arbitrary positive constants a, (3, 7, and c are related by a = P^c^', R = R/i? = 
(x — x')/|x — x'l; and the square brackets [ ] indicate that the enclosed quantity is to be 
evaluated at the retarded time f = t — R/c. 

To prove this theorem, we start by deriving the following identities from Eq. (1): 
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(4a) 
(4b) 



where c is arbitrary positive constant, 5 is the Dirac delta function, and p and J are arbitrary 
functions of space and time satisfying Eq. (1) at all points and at all times. 

Because Eq. (1) is satisfied at all points and at all times, we can evaluate it at the source 
point and the retarded time^ 



[V • J] + 



dp 
'dt 



0. 



(5) 
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A heuristic interpretation of Eq. (5) is as follows. Consider an observer at a particular 
location in space who has a watch that reads a particular time. The observer is surrounded 
by nested spheres, on each of which there is a well-defined retarded time (with respect to 
the observer). Equation (5) states that the continuity equation holds (or rather, held) on 
each of those spheres, at the relevant retarded time.^ 

We multiply Eq. (5) by and use d[p\/dt — [dp/dt] (see the Appendix) to obtain 



[V^ ■ J] 
R 



dt\R 



(6) 



We now substitute the identity (see the Appendix) 



[V^ • J] 

R 



into Eq. (6) and write 



[J 
R 



The gradient of Eq. (8) is 



V V 
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We now consider the identity (see the Appendix) : 



V V 



V X V X 



[J] 
R 



+ 



1 d^\J] 



Rc^ dt^ 

From Eqs. (9) and (10) and the property d[J]/dt = [dJ/dt], we obtain 
V X I V X 



-47r[J](5(x-x'). 



(9) 



(10) 



With the aid of the identities (see the Appendix) : 

V X 
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'dt 
dp 
di 
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we can see that Eq. (11) becomes Eq. (4a). To derive the identity (4b) we take the time 
derivative of Eq. (8) and use the property d[J]/dt — [dJ/dt] to obtain the expression 



We next substitute the identity (see the Appendix): 

1 ^2 
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(J^)+47rc^[p]5(x-x'). 



into Eq. (14) and obtain after rearranging terms 
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(16) 



Wc then substitute Eq. (13) into Eq. (16) to derive Eq. (4b). 

We next use Eq. (4) to obtain Eqs. (2) and (3). Note that the last term on the right-hand 
side in Eq. (4) , after integrated over all space, can be transformed into a surface integral that 
vanishes at infinity if the sources are localized, that is, when the sources are zero outside 
the surface of a finite region of space. Henceforth we will assume that p and J are localized 
sources.^ We multiply the second term on the left-hand side of Eq. (4a) by l3a/{Ana) and 
the remaining terms by /9/(47r) and integrate over all space: 
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The operators Vx and a /at can be extracted outside the integrals of the left-hand side. 
The first term on the right-hand side of Eq. (17) becomes AnPJ after integration over the 
delta function. The operator V can be extracted outside the last integral on the right-hand 
side, and the resulting integral can be transformed into a surface integral that vanishes at 
infinity because J is localized. Therefore, Eq. (17) reduces to 
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We next multiply Eq. (4b) by a /{An) and integrate over all space to obtain 
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If we follow an argument similar to that used to go from Eq. (4a) to Eq. (18), we obtain 



ap 
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(20) 



Equations (18) and (20) are the main result of this paper. These equations are actually 
Eqs. (3d) and (3a) because the quantities within the brackets { } in Eqs. (18) and (20) are 
the fields F and G defined by Eq. (2). In other words, Eqs. (18) and (20) show the existence 
of the fields F and G in terms of which these equations can be written as Eqs. (3d) and 
(3a). 

To complete the demonstration of the theorem we need to derive Eqs. (3b) and (3c). We 
first take the divergence of Eq. (2b) 



A direct calculation gives 
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(22a) 

(22b) 
(22c) 



where Eq. (AlO) and the results R • (A x R) = and V x (/(i?)R) = have been used. 
By a similar calculation we can show 



dJ_ 
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X 



R 



0. 



(23) 



If we substitute Eqs. (22) and (23) into Eq. (21), we obtain Eq. (3b). 
We now take the curl of Eq. (2a) 

dp 
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— f d'x'V X 



( 1 


'dJ' 




.9t_ 



A direct calculation gives 
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(24) 
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where Eq. (A12) and the results R x R = and V x (/(i?)R) = have been considered. 
By a similar calculation we can show 
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We now calculate 



V X 
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(27a) 
(27b) 

(27c) 



where Eq. (11) has been used. From Eqs. (24)-(27) and Eq. (2) we obtain Eq. (3c). This 
result completes the demonstration of the theorem. Note that another theorem that leads 
to Maxwell's equations has also been recently proposed. ^'^ 




III. MAXWELL'S EQUATIONS AND JEFIMENKO'S EQUATIONS 



If we identify p and J with the usual charge and current densities of electromagnetism 
and c with the speed of light in vacuum, then Eq. (2) become Jefimenko's equations^ 
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and Eq. (3) becomes Maxwell's equations, 



V • E = ap, 

V-B = 0, 

„ ^ dB 
VxE + ,- = 0, 

VxB-^— =/3J. 
a at 



(28a) 
(28b) 

(29a) 
(29b) 

(29c) 
(29d) 



Equations (28) and (29) are expressed in a generahzed system of units defined by a, /3, 
and 7 satisfying the relation 

a = f3-fc^. (30) 

This generalized system contains three systems of units: Gaussian, SI, and Heaviside- 
Lorentz. The specific values of a, /3 and 7 are given in Table I. If we adopt SI units, then 
Eqs. (28) and (29) become Jefimenko's equations^ and Maxwell's equations are expressed in 
their usual form.^ 
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System 


a 


/5 


7 


Gaussian 


in 


47r/c 


1/c 


SI 


1/eo 


/"o 


1 


Heaviside-Lorentz 


1 


1/c 


1/c 



TABLE I: Generalized system of electromagnetic units defined by a, l3, and 7 with a = /?7C^. This 
system contains the Gaussian, SI, and Heaviside-Lorentz units. 

IV. DISCUSSION 

The existence theorem for the retarded fields F and G in Eq. (2) is a purely mathematical 
result, that is, the theorem only ascribes mathematical existence to these fields by showing 
that they satisfy the field equations (3). To ascribe a physical meaning to these abstract 
fields, we must identify the conserved sources as the electric charge and current densities, 
and associate the fields with the electric and magnetic fields E and B, via the Lorentz force 
law, F = e(E + v/c x B). In other words, real electric and magnetic fields produce the 
observed forces on electric charges which are described by the Lorentz force. Therefore, the 
physical existence of the fields E and B can be verified by testing the Lorentz force, which 
itself constitutes an independent postulate. 

The theorem formulated here is formally correct. However, a referee pointed out that 
the physical implication of the theorem is not surprising and posed the following interesting 
question: what would it mean if the theorem were false from a physical point of view? The 
referee speculated that given the sources p and J satisfying the continuity equation, they 
would not yield the retarded fields E and B. We can imagine two situations. In one case we 
can conceive abstract charge and current densities p and J satisfying the continuity equation 
V ■ J + dp/ dt = which produce the instantaneous electric and magnetic fields (in Gaussian 
units) 

E(x,i)= J d'x'^p{^',t), (31a) 
B{^,t)=J dV:5Z^!^_^, (31b) 
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that satisfy the field equations of a Gahlean electromagnetic theory 



,12 



V-E 



47rp, 



(32a) 



V X E 



V B 



(32b) 
(32c) 




(32d) 



Our problem is that we do not know how to prepare the sources p and J in the laboratory. 
For example, we do not know how to produce sources p and J associated with a moving 
electron that would yield instantaneous fields. A moving (real) electron always produces 
retarded electric and magnetic fields satisfying Maxwell's equations, but not instantaneous 
electric and magnetic fields satisfying the Gahlean equations (32). We can also imagine 
that there exist physically realizable charge and current densities satisfying the continuity 
equation which do not generate electric and magnetic fields (at least, not fields satisfying 
Maxwell's equations). Because these situations are not possible, we conclude that given the 
conserved sources p and J, the fact that we can always construct electric and magnetic fields 
is not a surprise. 

V. CONCLUDING REMARKS 

The importance of charge conservation was noted by Maxwell who used it in his discovery 
of the displacement current. This current allowed him to extend Ampere's law (formulated 
for a steady-state regime) to the non- stationary regime. Since then, most textbooks intro- 
duce the displacement current by invoking charge conservation. But if charge conservation 
is used to find the final form of Maxwell's equations, then the alert reader might find it 
surprising that this conservation law can also be considered as a consequence of Maxwell's 
equations — the standard proof being that Maxwell's equations directly imply the continuity 
equation. The idea that charge conservation is not an independent assumption, but a con- 
sequence of the laws of electrodynamics is recurrent in textbooks. The reader can see here 
a circular argument: the continuity equation can be obtained from Maxwell's equations but 
in obtaining Maxwell's equations, the continuity equation is usually assumed. Our point of 
view is that charge conservation can be considered as a postulate rather than a result of the 
theory. Our formal approach to Maxwell's equations based on the existence theorem shows 
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an example of how a single postulate (continuity equation) of a field theory may be used to 
derive not only the retarded fields associated with that theory, but also the explicit form of 
the field equations satisfying such fields. 
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APPENDIX: IDENTITIES INVOLVING RETARDED QUANTITIES 

We present a simple and non-rigorous derivation of some identities for retarded functions. 
A retarded function is denoted by [p] and is defined as [p] = p{x.',t — R/c), where c is an 
arbitrary positive constant, or equivalently, by [p] = p{x.',t'), where t' — t — R/c. We note 
that [p] depends on the source coordinates not only explicitly, but also implicitly through 
R. Therefore, the derivation of identities involving space derivatives of retarded quantities 
is a somewhat complicated task. 

We can represent a retarded quantity using the Dirac delta function: 



where u — t' — t + R/c. This representation of a retarded quantity applies also to vector 



for example, [V'p] = Jdt' S{u)V' p{^' ,t'); [dp/dt] = Jdt' 6{u)dp{x',t')/dt' and [V • J] = 
Jdt'S{u)V' ■ J{x',t'). We will now obtain Eqs. (7), (10), (12), (13), and (15). We will 




(A.l) 



functions: [J] = / dt' 5(xi) J(x', t'). Similar representations can be written for other functions. 



require the following derivatives involving the delta function: 



dSju) 
dt 



d5{u) 

-V'5(u), 
R dSju) 

1 d^ju) _ 2 dSju) 

c2~at2 r2 dt 



(A.2) 



(A.3) 



(A.4) 



(A.5) 
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We begin by taking the derivative of [p] with respect to the present time: 



dt 





J 


1 dtS{u) 




dp 


-6{u)p{^',t') 


t'=-oo 
t'=-oo 


dp 






at 







(A.6a) 
(A.6b) 



(A.6c) 

where we have used Eq. (A2) and done an integration by parts using the result 5(±oo) = 0. 
By a similar calculation we can show d[i\/dt — [di/dt\. 

The divergence of [J] with respect to the source coordinates gives 



' • [J] dt'V'5{u) ■ J(x',0 + j dt'5(u)V ■ J(x',i') 
= -J dt'VS{u) ■ J(x',i') + J dt'6{u)V ■ J(x',i' 
= -V • [J] + [V ■ J] , 



where we have used Eq. (A3). By combining the expansions: 

v.H4v.[Jl + [Jl.v(i 



v'.l^4v'.iJi + iJi.v'(i). 



(A.7a) 
(A.7b) 
(A.7c) 

(A.8) 
(A.9) 



and using Eq. (A7) we obtain Eq. (7). 

The curl of [J] with respect to the field coordinates is 



V X fJl 



J dt'V6{u) X J(x',t') 
j dt'5{ 



_R d 



M J X',f 



' c"" at 



d3_ 
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X 



(A.lOa) 
(A.lOb) 
(A. 10c) 



where Eq. (A4) and d[J]/dt = [dJ/dt] have been used. We use Eq. (AlO) and obtain Eq. (12) 



V X 
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iv X [J] - [J] X V (i) 
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The gradient of [p] with respect to the field coordinates is 

V[p] = 



dt'V6{u)p{x',t') 

dt'S{u)p{^',t') 



R d 



R 

c 



dp 
di 



where Eqs. (A4) and (A6) have been used. From Eq. (A12) it follows that 
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The Laplacian of [J] with respect to the field coordinates is 

v'[j] =y" dtv^s{u)j{K',t') 



1^ 

1 d^[J] 2 d[J] 



[ dt'5(u)J(:>c',t') - [ dt'5(u)J(K',t') 

J RcdtJ 



(A.12a) 
(A. 12b) 
(A. 12c) 

(A.13) 

(A.14a) 
(A. 14b) 
(A.14c) 



c2 dt^ Rc dt ' 

where Eq. (A5) has been used. Evidently, Eq. (A14) is vahd for the components of J. For 
example, if Jx is the component x of J, then 



1 dVx] 2 d[Jx 
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Rc dt 



(A.15) 



Consider now the Laplacian of [Jx\/R with respect to the field coordinates: 
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1 d^[Jx 
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R 



R^c dt 
-47r[J,](5(x-x'), 



fid[Jx 
c dt 



R 



(A.16a) 

- 47r[Ja;]5(x - x') (A.16b) 
(A. 16c) 



where Eq. (A12) [with instead of p] and Eq. (A15) have been used. Similar expressions 
for the components Jy/R and Jz/R can be found, and therefore we can write the identity 
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1 d^\J] 



R J Rc^ dt^ 
which yields Eq. (10) after using the expansion 

Eq. (A16) [with p instead of Jx] directly implies Eq. (15). 



= -47r[J](5(x-x'), (A.17) 
= V(V • ) - V X (V X ). Finally, 
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